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Near comoving wavenumber k, the gravitational-wave background (GWB) from inflation carries 
information about the physical conditions near two moments in cosmic history: the moment when k 
“left the horizon” during inflation, and the moment when it “re-entered the horizon” after inflation. 
We investigate the extent to which this information can be extracted if the GWB is measured by a 
combination of cosmic-microwave-background (CMB) polarization experiments on large scales and 
space-based laser-interferometer experiments on small scales. To disentangle this information, we 
derive a new gravitational-wave transfer function that incorporates a number of physical effects 
that were treated less accurately, less generally, or were missing altogether in previous treatments. 
In particular, it incorporates: (i) dark energy with time-varying equation-of-state w{z); (ii) tensor 
anisotropic stress due to free-streaming relativistic particles in the early universe; and (Hi) a variety 
of physical effects that cause deviations from the standard equation-of-state w = 1/3 during the 
radiation era. Based on this transfer function, we consider the degree to which the GWB can be 
used to test inflation and to probe the “primordial dark age” between the end of inflation and the 
electroweak phase transition. 


I. INTRODUCTION 

Inflation P| generates tensor perturbations (gravita¬ 
tional waves) with a nearly scale-invariant primordial 
power spectrum . This gravitational-wave background 
(GWB) contributes to the cosmic microwave background 
(CMB) temperature anisotropy ill, and also produces 
a characteristic curl ^r “B-mode”) fingerprint in the 
CMB polarization 0, IE 0 ■ Several CMB experiments 
are being developed to pursue this signal Q- The GWB 
also persists as a sea of re lic g ravitational radiation fill¬ 
ing the universe today ■ Direct detection of this 

relic radiation has received considerable attention over 
the past year or so, since it has been realized that space- 
based laser interferometers operating in the frequency 
range 0.1 Hz < / < 1 Hz might achieve the necessary 
sensitivity and foreground subtraction pm In par¬ 
ticular, two futuristic experiments have been proposed — 
NASA’s “Big Bang Observer” (BBO) and the Japanese 
“Deci-hertz Interferometer Gravitational Wave Observa- 
tory” (DECIGO) — and are currently being investigated 

PPP. 

The gravitational wave spectrum generated by infla¬ 
tion carries important information about the conditions 
during inflation. But the spectrum also receives correc¬ 
tions, both large and small, from the subsequent evolu¬ 
tion and matter content of the universe after inflation. 
In this paper, we identify various post-inflationary phys¬ 
ical effects which modify the GWB, and show how they 
may be encoded in the gravitational-wave transfer func¬ 
tion that relates the primordial tensor power spectrum to 
the gravitational-wave spectrum at a later point in cos¬ 
mic history. It is necessary to properly understand and 
disentangle the post-inflationary effects in order to opti¬ 
mally extract the inflationary information in the GWB. 
But these modifications are also interesting in their own 


right, since they offer a rare window onto the physical 
properties of the high-energy universe during the “pri¬ 
mordial dark age” between the end of inflation and the 
electroweak phase transition. 

The same features that make the inflationary GWB 
difficult to detect — namely its small amplitude and 
gravitational-strength coupling to matter — also make 
it a clean cosmological probe. First, because of their tiny 
amplitude, the gravitational waves obey linear equations 
of motion, so that their evolution may be predicted an¬ 
alytically with high precision. (By contrast, density per¬ 
turbations grow after horizon entry, and perturbations 
shorter than ^ 10 Mpc have already grown non-linear; 
so it is impossible to study their evolution analytically, 
and even numerical analysis is challenging.) Second, be¬ 
cause of their ultra-weak interactions with matter, the 
gravitational waves have been free streaming since the 
end of inflation — in contrast to neutrinos (which began 
streaming roughly a second later) and photons (which 
began streaming several hundred thousand years later). 
The gravitational waves carry unsullied information from 
the early universe, and subsequent modifications of the 
gravitational-wave spectrum are not washed out by ther¬ 
mal effects (since the gravitons are thermally decoupled). 

The gravitational-wave spectrum near a given 
wavenumber k is primarily sensitive to two “moments” 
in cosmic history: (I) the moment when that mode “left 
the horizon” (i.e., became longer than the instantaneous 
Hubble radius during inflation), and (2) the moment 
when the mode “re-entered the horizon” {i.e., became 
shorter than the instantaneous Hubble radius once again, 
after the end of inflation). The first moment imprints in¬ 
formation about inflation itself, while the second moment 
imprints information about post-inflationary conditions. 
The GMB is sensitive to long-wavelength modes that re¬ 
entered at relatively low temperatures (well after BBN), 
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corresponding to relatively well-understood physics. By 
contrast, laser interferometers are sensitive to shorter 
wavelengths that entered the horizon at high tempera¬ 
tures (T ~ lO”^ GeV), well above the electroweak phase 
transition. The physical conditions at such high energies, 
which are considerably beyond the reach of particle ac¬ 
celerators, are a mystery, so that any information about 
this epoch from the GWB would be very valuable. 


The outline of our paper is as follows. In section nn 
we review the basics of inflationary gravitational waves. 
We take special care to clearly establish our conventions, 
and to point out where conventions diverge and become 
confused in the literature. In section IIIII we present a 
new calculation of the gravitational-wave transfer func¬ 
tion which includes effects not considered in previous cal¬ 
culations [Hill HI We separate the transfer function 
into three factors, each with a distinct physical meaning. 


The first factor accounts for the redshift-suppression of 
the gravitational-wave amplitude after horizon re-entry. 
Among other things, this factor accomodates a dark en¬ 
ergy component with a time-varying equation-of-state 
parameter w{z). The second factor captures the be¬ 
havior of the background equation-of-state parameter w 
near the time of horizon re-entry. The third factor ac¬ 
counts for the damping of tensor modes due to tensor 
anisotropic stress from free-streaming relativistic parti¬ 
cles in the early universe. This damping effect was re¬ 
cently pointed out by Weinberg [a, and the damping 
on CMB scales due to free-streaming neutrinos has been 
studied in [Hlllllllll In this paper, we stress that 
it is also necessary to consider this damping effect on 
laser-interferometer scales, which re-entered the horizon 
when free-streaming particles were an unknown fraction 
/ of the background energy density. We present accurate 
expressions for the damping effect as a function of /. We 
also observe that Weinberg’s analysis, which originally 
focused on a single fermionic particle (the neutrino), ex¬ 
tends in a simple way to the more general case of a mix¬ 
ture of free-streaming bosons and fermions with different 
temperatures and decoupling times. 


We identify six physical effects which can modify the 
relic GWB by causing the equation-of-state parameter 
w to deviate from its standard value {w = 1/3) during 
the radiation-dominated epoch. Furthermore, although 
it is often treated as a stationary random process, the 
inflationary GWB is actually highly non-stationary (as 
emphasized by Grishchuk j^; also see IHES)- Thus, 
our transfer function keeps track of the coherent phase 
information that it contains. Finally, in section CYI we 
discuss the implications of these results for future CMB 
and laser-interferometer experiments, and consider what 
we might learn by measuring the GWB in both ways. 


II. TENSOR PERTURBATIONS: 
FUNDAMENTALS AND CONVENTIONS 


In this section, we derive some basic facts about in¬ 
flationary gravitational waves. These derivations are not 
new, but are intended to make the present paper more 
pedagogical and self-contained. They also establish our 
conventions explicitly, and provide a brief guide to con¬ 
ventions used by other authors. Such a guide is neces¬ 
sary since there are many slightly different conventions 
floating around in the inflationary gravitational-wave lit¬ 
erature and, as a result, erroneous factors of 2 and tt are 
ubiquitous. We have been careful to highlight each spot 
where convention choices crop up, and to state which 
convention we have chosen. 

Tensor perturbations in a spatially-flat Friedmann- 
Robertson-Walker (FRW) universe are described by the 
line element 


ds^ = a^[—dr^ -|- {5ij + /ly )da;Ma:'^], (1) 

where t is the conformal time, cc® are comoving spatial 
coordinates, and /ly is the gauge-invariant tensor metric 
perturbation. Note that, although our definition of /ly is 
perhaps the most common (see e.g. some au¬ 

thors define the tensor perturbation with an extra factor 
of 2 (see e.g. |HI13) that the spatial metric is written 
Sij + 2hij. Throughout this appendix, if a term contains 
a repeated index, that index should be summed: from 1 
to 3 for latin indices and from 0 to 3 for greek indices. 
The perturbation /ly is symmetric (/ly = hji), traceless 
{hii = 0), and transverse {hijj = 0), and therefore con¬ 
tains 6 — I — 3 = 2 independent modes (corresponding to 
the “-I-” and “x” gravitational-wave polarizations). 

One can think of /ly (r, x) as a quantum field 
in an unperturbed FRW background metric = 
dfagi—a^, a^, a^, a^}. At quadratic order in /ly (which 
is adequate, since /ly is tiny), tensor perturbations are 
governed by the action (see e.g. jHIlll) 


S = 


dT(hiy/—g 


' 1 
647rG ^ 2 ^ ^ 


( 2 ) 


where g^'^ and g are the inverse and determinant of 
respectively, and G is Newton’s constant. The tensor 
part of the anisotropic stress 

Ily = Tj — p6^j (jp = unperturbed pressure) (3) 

satisfies 11^^ = 0 and djlij = 0, and couples to hij like an 
external source in ©• By varying hij in ©, we obtain 
the equation of motion 

h'h + 2^/iT - V^/iy = 167rGa2(r)ny (r, x), (4) 

where a prime (') indicates a conformal time derivative 
d/dr. Next, it is convenient to Fourier transform as fol- 
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lows: 


KjiT.yi) = ^ V IQ-kG I^^^ 3/2 (k)fek(T)e^'^^, (5a) 

n,^(r,x) =^Vl^Jj^^el^{k)UUr)e^-, (5b) 

where r =(“+” or “x”) labels the polarization state, and 
the polarization tensors are symmetric [e(),(k) = ejj(k)], 
traceless [eii(k) = 0], and transverse [kie\j{k) = 0]. We 
also choose a circular-polarization basis in which (k) = 
{e^j (-k))* , and normalize the polarization basis as follows: 


^e-.(k)(e^^.(k))*=2J-. (6) 


Although our normalization convention is the most 
standard one, other conventions — i.e. different numeri¬ 
cal constants on the right-hand side of — exist in the 
literature. Substituting m into then yields 


S='^fdTdk^ 


hl'h!' - /if -H327rGa^nr/r 






( 7 ) 


Now we can canonically quantize by promoting /ij) and 
its conjugate momentum 

^l{T)=a\rW^{r) (8) 


to operators, /ij) and 7rj), satisfying the equal-time com¬ 
mutation relations 


= i^-<5(3)(k_k'), (9a) 

Kir),hl{r)] = K(r),7r£,(r)]=0. (9b) 


Since hij (r, x) is hermitian, its fourier components satisfy 
~ write them as 

= hf,{T)al + hl{T)a^^, ( 10 ) 


where the creation and annihilation operators, and 
aj), satisfy standard commutation relations 





(11a) 






= 0 , 


(11b) 


while the (c-number) mode functions hk (t) and /i^ (r) are 
linearly-independent solutions of the fourier-transformed 
equation of motion 

K + 'i'^K + t^h, = W,GaHr)UM. ( 12 ) 

Eq. (irnil makes use of the fact that, by isotropy, the 
mode functions hf,{T) will depend on the time (r) and the 


wavenumber (k = |k|), but not on the direction (k) or the 
polarization (r). Note that consistency between the two 
sets of commutation relations, m and CD, requires that 
the mode functions satisfy the Wronskian normalization 
condition 


Ki^Wkij) - hl(T)h'^{T) = (13) 

a?{T) 

in the past. In particular, the standard initial condition 
for the mode function in the far past (when the mode k 
was still far inside the horizon during inflation): 




exp(—i/cr) 
a(r)-\/2fc 


(as r ^ —oo), (14) 


satisfies — but it is not the unique initial condi¬ 
tion which does so. This is a manifestation of the well 
known vacuum ambiguity that is responsible for particle 
production in cosmological spacetimes (see 

Now that we have discussed the quantization of tensor 
perturbations, let us turn to the three different spectra 
that are commonly used to describe the stochastic GWB: 
the tensor power spectrum A^(/c, r), the chirp amplitude 
/i^(fc,r), and the energy spectrum flg^(k,T). 

In the early universe, the GWB is usually character¬ 
ized by the tensor power spectrum A^(fc,r). With the 
formalism developed thus far, one can check that 


r°° k^ ^dk 

(0|/iy(T,x)/i,j(T,x)|0)=y^ 647rG—|/ife(T)| (15) 

so that the tensor power spectrum is given by 


Aft(fc,r) 


d{0\hjg\0) 

dink 


= 647rG^ \hk{T)\^ ■ 


(16) 


Note that our definition fTC|) agrees with the WMAP con¬ 
vention (see [29)1 — this will be clearer when we present 
the approximate slow-roll form of the spectrum below. 
Although the WMAP convention seems to be becoming 
the standard one, several other definitions of the tensor 
power spectrum exist in the literature, and differ from 
(Cnj by an overall numerical constant. Also, since itTO) 
defines the tensor power spectrum in terms of the full 
tensor perturbation hij^ the normalization of the power 
spectrum is independent of the normalization of the 
polarization basis. By contrast, some authors define the 
tensor power spectrum in terms of the polarization com¬ 
ponents of hij , so that the normalization of the spectrum 
is linked to the convention-dependent coefficient on the 
right-hand side of ®- 

The present-day GWB is usually characterized either 
by its chirp amplitude h^{k,T), or by its energy spec¬ 
trum flg^{k,T). First, the chirp amplitude represents 
the rms dimensionless strain (^ dL/L in a gravitational 
wave antenna) per logarithmic wavenumber interval (or 
logarithmic frequency interval), and is related to the ten¬ 
sor power spectrum as follows: /i^(fc,r) = A/A((fc,r)/2 
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(see 1^ for more details). The energy spectrum 

1 


^qw — 


PcritiT) d\nk 


(17) 


represents the gravitational-wave energy density (pg^) 
per logarithmic wavenumber interval, in units of the 
“critical density” 


Peril ('^) 


3g^(r) 

SttG 


(18) 


To compute ilg^{k,T), we can again think of the tensor 
perturbation hij as a field in an unperturbed FRW back¬ 
ground metric and then use its action |(2Jl to compute 
its stress-energy tensor 

Taf3 = —+ PapL, (19) 

where L is the Lag rang ian function in @ — see, for ex¬ 
ample, Ch. 21.3 in |31l . In particular, the gravitational- 
wave energy density (ignoring the anisotropic stress cou¬ 
pling) is 


where, in this equation (and in the remainder of this 
section) “=” denotes a definition, indicates that 

the slow-roll approximation has been used, and Mpi = 
(87rG)“^/^ is the “reduced Planck mass.” Note that the 
primordial power spectrum is time-independent, since 
(by definition) it is evaluated when the mode k is far out¬ 
side the horizon. Although our derivation of may 
seem crude, it is well known that (1^ provides a very 
good approximation to the inflationary tensor spectrum. 
We will not reproduce it here, but a closely analogous 
derivation leads to the slow-roll expression for the pri¬ 
mordial scalar power spectrum: 


A^(fc) 


1 




2ttM, 


pi 


2 


(25) 


It is also useful to define the tensor/scalar ratio 


r{k) = Al{k)/A^{k) = 16e,, (26) 


along with the scalar and tensor spectral indices 


1 ihP)^ + iWh.gr 

° 647rG a2 

which has vacuum expectation value 

k^ + k"^ dk 


(0|p,„|0)= [ 

Jo 


2tt^ o2 k ’ 

so that the energy spectrum is given by 
8ttG k^ 


^qwik,T) = 


3772 (r) 27r2 


rr) 


( 20 ) 


( 21 ) 


( 22 ) 


We have shown how each of the three spectra (Aj, 
and ) may be written in terms of the mode functions 
/ij,(r). Also note that, once the mode k re-enters the 
horizon after inflation, the corresponding mode function 
obeys = fc^|/ifc(T)p, as explained in the next 

section, so that the three spectra may be related to one 
another in a simple way: 


1 fc^A2(fc,T) 1 fc^/i^(fc,r) 

12 a2(T)772(7-) 6 a 2 (T) 772 ( 7 -( 


(23) 


Next, let us “derive” the slow-roll expression for the 
primordial tensor power spectrum. As long as k remains 
inside the Hubble horizon (k ^ aH) during inflation, 
the mode function (t) is given by (II3; and once k is 
outside the horizon {k <C aH), the mode function 
is independent of t. Then, by simply matching 
to at the moment of horizon-exit {k = aH), one 

finds = l/(a^\/^), where an asterisk (*) denotes 

that a quantity is evaluated at horizon-exit (fc = a*77^). 
Thus, the primordial tensor power spectrum is given by 


/ N dlnA^ 

n.(k) . 


d\nk 


-6e* - 

-2e.,. 


2??*, 


(27a) 

(27b) 


In the slow-roll formulae above, we have introduced the 
usual “potential” slow roll parameters 




= \m. 


,-..2 V'i^) 


r] = M 


pi' 


v{4>) 


(28) 


where V{(/)) is the inflaton potential. Note that Eqs. 12611 
and (|27b|) together imply the well known consistency re¬ 
lation for inflation 


r = —Srif.. (29) 

III. THE TENSOR TRANSFER FUNCTION 

Since cosmological tensor perturbations are tiny, they 
are well described by linear perturbation theory, so that 
each fourier mode evolves independently. Thus, we see 
from Eq. 03 that the primordial tensor power spectrum 
— defined at some conformal time shortly after the 
end of inflation, when all modes of interest have already 
left the horizon, but have not yet re-entered — is related 
to the tensor power spectrum at a later time r by a mul¬ 
tiplicative transfer function 

Kik, t) = Tf.{k, T)Al{k, Ti), (30) 


where 


Aiik) = 647rG 


27r2 


Tout) 


77* 


(^2Tr Mpi 


(24) 


Th{k,T) 


hkjr) 

hk{T.) 


( 31 ) 
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Note that we will not necessarily want to evaluate 
Th(k,T) at the present time (r = Tq), since different ex¬ 
periments probe the gravitational-wave spectrum at dif¬ 
ferent redshifts. For example, while laser interferometers 
measure today, CMB experiments measure it near 
the redshift of recombination. As long as a mode re¬ 
mains outside the horizon (k<^aH), the corresponding 
perturbation does not vary with time [h^.{T) = hf^{Ti)\, so 
that the transfer function is very well approximated by 
Th{k,T) = 1. (For a general proof, even in the presence 
of anisotropic stress, see the appendix in M-) Thus, the 
rest of this paper will focus on T^{k^T) for modes that 
have already re-entered the horizon prior to time r. 

It is very convenient to split the transfer function m 
into three factors as follows: 


Th{k,T) 


hkij) hkij) h^jr) 
hk{T^)hk{T) hk{T) 


C1C2C3. 


(32) 


Here hkir), /ife(r) and hkir) represent three different so¬ 
lutions of the tensor mode equation m- In particu¬ 
lar, /ifc(r) is the true (exact) solution of HI211 : hkir) is 
an approximate solution obtained by ignoring the ten¬ 
sor anisotropic stress Hfc on the right-hand-side of m-, 
and /ifc(r) is an even cruder approximation obtained by 
first ignoring Hfc and then using the “thin-horizon” ap¬ 
proximation, described in subsection IIII Al to solve (I12II . 
[Briefly, the thin-horizon approximation treats horizon 
re-entry as a “sudden” or instantaneous event. In this 
approximation, hki^r) is frozen outside the Hubble hori¬ 
zon, redshifts as l/a(T) inside the Hubble horizon, and 
has a sharp transition between these two behaviors at 
the moment when the mode re-enters the Hubble hori¬ 
zon (fc = aH).] 

These three factors each represent a different physical 
effect. The first factor, 


Cl 


hkir) 

hk{Ti) 


(33) 


full transfer function Th- The other two factors, C 2 and 
C 3 , are typically of order unity, and are primarily sen¬ 
sitive to the physical conditions near the time that the 
mode k re-entered the Hubble horizon. In subsections 
iTTnniTTm and irrm we derive expressions for Ci, C 2 
and C 3 , respectively, and flesh out the physical inter¬ 
pretations given above. In subsection HTTPl we discuss 
various physical effects that cause deviations 5w from 
the usual equation-of-state parameter {w = 1/3) in the 
early universe, and explain how these effects modify the 
gravitational-wave transfer function. 


A. The redshift-suppression factor, Ci 

The mode function h^(T) behaves simply in two 
regimes: far outside the horizon (fc ^ aH), and far 
inside the horizon {k ^ aH). Far outside, /ij,(t) is r- 
independent (as we have seen). Far inside, after horizon 
re-entry, hf.{T) oscillates with a decaying envelope 

= ^^cos[A:(T-rfc)-t(/fc], (36) 

as we shall see in the next two subsections, where A^, and 
are constants representing the amplitude and phase 
shift of the oscillation, and is the conformal time at 
horizon re-entry (k = aH). These two simple regimes are 
separated by an intermediate period (horizon-crossing) 
when k ^ aH. 

In the thin-horizon approximation, we neglect this in¬ 
termediate regime. That is, we assume that hj^ir) = 
hkir^) when k < aH-, and that h^(r) is given by Eq. 
(TTHli for k > aH ; and that the outside amplitude is con¬ 
nected to the inside envelope via the matching condition 

('’’*) = Ak/a{Tk). Ignoring the phase shift (pk, which 
is really an asymptotic quantity, this matching condition 
simply imposes continuity of the inside and outside am¬ 
plitudes at k = aH. Combining the matching condition 
with Eq. (ESI, we see that 


represents the redshift-suppression of the gravitational- 
wave amplitude after the mode k re-enters the horizon. 
The second factor. 


C 2 


hk{T) 

hkir) 


(34) 


captures the behavior of the background equation-of- 
state parameter w{t) = p(r)/p(r) around the time of 
horizon re-entry. The third factor. 


C 3 


hkjr) 

hkir) 


(35) 


measures the damping of the gravitational-wave spec¬ 
trum due to tensor anisotropic stress. Note that Ci by 
itself is <C 1 and provides a rough approximation to the 


Cl = cos^[k(T - Tk) -b (/fc], (37) 

where 1 + z = aQ/o(T) is the redshift at which the spec¬ 
trum is to be probed, and 1 +^^ = a^/a^ is the redshift at 
which the mode re-entered the Hubble horizon (k = aH). 

The relic GWB from inflation is often treated as a 
“quasi-stationary” process (which means that its statis¬ 
tical properties only vary on cosmological time scales — 
much longer than the timescales in a terrestrial experi¬ 
ment). But the cos^[...] factor in Eq. (I37II implies that 
the background is actually highly non-stationary — its 
power spectrum oscillates as a function of both wavenum¬ 
ber k and time t. This cos^[...] factor represents a gen¬ 
uine feature, and is not a spurious byproduct of our thin- 
horizon approximation. Physically (as observed in 
the inflationary GWB consists of gravitational standing 
waves with random spatial phases, and coherent temporal 
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phases. All modes k at fixed wavenumber k = |A:| re-enter 
the Hubble horizon simultaneously, and subsequently os¬ 
cillate in phase with one another — even until the present 
day. Classically, the modes are synchronized by inflation; 
quantum mechanically, they are squeezed |^. Thus, at 
a fixed wavenumber, Th{k,T) is sinusoidal in r, with os¬ 
cillation frequency A:, and a phase shift (fij, (computed 
in the next subsection). Alternatively, at fixed time, 
Th{k,T) oscillates rapidly in wavenumber. For the direct 
detection experiments we are considering, which cannot 
resolve these oscillations, the factor cos^[...] should be 
averaged, and replaced by 1/2 in Eq. iHTTll . 

In the remainder of this subsection, we derive an accu¬ 
rate expression for (l-|-z^), the horizon-crossing redshift. 
To start, let us write the background energy density p as 
a sum of several components. The Ah component has en¬ 
ergy density Pj, pressure equation-of-state parameter 
Wj = p ^/and obeys a continuity equation 

^=3[l + zn,(.)]^. (38) 

Pi 1 i" 2 


Integrating this equation yields 




(1 -f z)^exp 



w,{z) 
1 + 2 



(39) 


where pf''^ is the present value. Then the Friedmann 
equation 




(40) 


may be rewritten as 


2jj2 X! ^ 


( 0 ) 


(1 -I- 2 :)exp 


Wjjz) 

1 -f 5' 


dz 


(41) 


where Hq is the present Hubble rate, and the density pa¬ 
rameter = p^P/p^cr represents the Ah component’s 
present energy density in units of the present critical den¬ 
sity p^cr = SHq/8ttG. Hence z^. is obtained by solving 
the equation 


(k/kof = F{z,) 


where 
Fi^k) = 

i 



w,{z) 
1 + z 



(42) 


(43) 


and fcp = QqHq = h x 3.24 x 10“^® Hz is today’s horizon 
wavenumber. 

Before solving this equation properly, let us pause to 
extract a few familiar approximate scalings from our for¬ 
malism. Since the primordial inflationary power spec¬ 
trum A^(fc,Ti) is roughly scale invariant [cx (A:/A:q)°], the 
current power spectrum A^(fc,rg) is roughly cx Ci, and 
hence cx (1 -I- Zf.)~'^. From Eq. (EHll we have h^ik^T^) cx 


(1 + ^fc) ^ and flg^ik^To) cx {k/ko)'^{l + z^.) For 
modes that re-enter the horizon during radiation dom¬ 
ination, when the Wj. = 1/3 term dominates the sum 
m, we solve to find (1 -I- Zf.) oc (fc/fco), which im¬ 
plies the approximate scalings h^{k,To) oc (k/ko)~^ and 
^gwik,To) cx (k/ko)^. For modes that re-enter during 
matter domination, when the = 0 term dominates 
the sum we find (1 + Zk) cx (fc/feg)^, which implies 
that hJJt,TQ) and flg^{k,To) are both cx (fc/fcg)“^ in this 
regime. 

For a more proper analysis, consider a universe with 4 
components: matter = 0), curvature (w^ = —1/3), 
dark energy (?w^g(z)), and radiation {w^{z) = 1/3 -I- 
dWj.{z)). Note that, although one often assumes = 1/3 
during radiation domination, we have allowed for cor¬ 
rections Sw^{z) due to early-universe effects discussed in 
subsection inmi Then we can write 

F{z) = F{z) + 6F{z), (44) 


where 


F{z) = 0/°) (1 + zf + (1 + z) + np , (45) 


and 


SF{z) = H|i°j(l-hz)exp 




HW(l+z)2{exp 


l-k5 J 
Sw^{z) 


3 dz 

Jq 1 + ^ 


-1 


}. (46) 


Here F represents a universe with spatial curvature, mat¬ 
ter, and “standard” {w^ = 1/3) radiation; and SF con¬ 
tains the modifications due to dark energy (w^^) and 
equation-of-state corrections (Sw^). 

If we neglect these modifications [by setting 12^ = 0 = 
Swr(z) so that SF = 0], Eq. B^ has the exact solution 


1 + 2/^ = 


“eg 


-1 




1 -k 


4i(k/ko)^ - nP 


(1 + ^eg)^' 


K 

( 0 ) 


(47) 


where 1 -I- z^g = fiP/flP is the redshift of matter- 
radiation equality. Then, including both modifications, 
the solution becomes 


(1 + ^fe) — (1 + Zk) + 

where is defined by Eq. EH) and 




F'(Sk) \ . I /, ^ F"(z,)SF(z'J 

F"(2,) I ^ V lF'(z,)P 


(48) 


(49) 


This solution is obtained by Taylor expanding F{zf/) 
around Zf. (to 2nd order in Szj. = Zf^ — z^.), and then 
solving the equation E(z^) = F{zf/) for 5z^. It is ex¬ 
tremely accurate for a wide range of w^^{z), and 
Sw^{z). Indeed, the simpler Ist-order expression 


5F{2k) 


F'{2^) 


( 50 ) 


is often sufficiently accurate. 
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B. The horizon-crossing factor, C2 

In the previous subsection, we treated horizon re-entry 
as a sudden event that occurs when k = aH. In reality, 
the “outside” behavior {hk = constant) only holds when 
k aH, and the “inside” behavior {hk oc a“^cos[fcr -|- 
phase]) only holds when k ^ aH. In between, when 
k ~ aH, neither behavior holds — i.e., the horizon has 
a non-zero “thickness.” 

The behavior of the background equation-of-state 
w{t) = p{t)/p{t) during the period of horizon re-entry 
is imprinted in the factor 02- For example, let us com¬ 
pute C 2 for a mode k that re-enters the Hubble horizon 
when w{t) is varying slowly relative to the instantaneous 
Hubble rate. Then we can write 

2 

a = ao(r/ro)“ with a = ^ , (51) 

so that the equation of motion for hk 

h'l+ 2{a'/a)h'i. +k’^h). = 0 (52) 

has the solution 

hk{T) = hk{T,)T{a + l/2)[Ar/2]i/2-«4_^^2(fcr), (53) 

where we have used fcrj<Cl and /i5,(rj) = 0. (Early in the 
radiation era, the relevant modes were far outside the 
horizon, and the corresponding mode functions were r- 
independent.) We have neglected the spatial curvature, 
K, because the two conditions K a^H^ (current ob¬ 
servations indicate that the spatial curvature is small) 
and k > aoHo (we are only interested in modes that are 
already inside the horizon) imply that K produces a neg¬ 
ligible correction to the equation of motion for hk- 
Once the modes are well inside the horizon {kr ^ I), 
we can use the asymptotic Bessel formula 

Ja-i/ 2 {kT) cos(fcT - a7r/2) (54) 

to find 

^ _ F ^^.^^ 2 ]~^'^cos^(fcr — ( 55 ) 

On the other hand, since a mode k re-enters the horizon 
{k = aH = a'/a) at time Tk = ajk, we can rewrite Eq. 
(TTtIi for C\ as 

Cl = [T/rfc]^“cos^[fc(T - Tk) + (j)k], (56a) 

= [/cr/a]“^“cos^(fcr — a -I- <j)k). (56b) 

Comparing Eqs. ( 1 ^ . (IMIl . and Il56bll . we see that 
the phase shift (j)k Ed- m is given by 

4 = [1 - 7r/2]Q;, (57) 


Mi 

Ki' 


where a should be evaluated at horizon re-entry (fc = 
aH). In particular, note that 

w=0 C 2 {k) = ^ and <1)^ = 2-tt, , . 

w = ^ => C 2 {k) = 1 and t/fj. = I — tt/2. ^ 


C. The anisotropic-stress damping factor, C3 


In this subsection, we will include the effects of the 
anisotropic stress term H^ on the right-hand side of 
the tensor mode equation m- A non-negligible tensor 
anisotropic stress H^ is most naturally generated by rel¬ 
ativistic particles free-streaming along geodesics that are 
perturbed by the presence of tensor metric perturbations 
hj,. In this situation, Weinberg has recently shown 
that the tensor mode equation 03 may be rewritten as 
a fairly sim ple integro-differential equation for hf. — see 
Eq. (18) in [ 13 . 

Let us focus on a particularly interesting case: a 
radiation-dominated universe in which the free-streaming 
particles constitute a nearly-constant fraction / of the 
background (critical) energy density. (Physically, if the 
free-streaming particles are stable, or long-lived relative 
to the instantaneous Hubble time at re-entry, then / will 
indeed be nearly-constant, as required.) In this case, fol¬ 
lowing an approach that is essentially identical to the one 
outlined in jlil . we write the solution in the form 

00 

hkir) = hk{T,) aj^ikr), (60) 

n—0 

where j^{kT) are spherical Bessel functions, and find the 
first five non-vanishing coefficients to be given by 


10 / 

"(15+ 4/)’ 

18/(3/+ 5) 

(15 + 4/)(50 + 4/)’ 

if2/(14/2 + 92/ + 35) 

(15 + 4/)(50 + 4/)(105 + 4/)’ 
^/(4802/3-H78266/2-hl61525/-29400) 
^ (15 + 4/)(50 + 4/)(105 + 4/)(180 + 4/) 


(61a) 

(61b) 

(61c) 

(61d) 

(61e) 


The odd coefficients all vanish: a 2 n+i = 0. Keeping these 
first five non-vanishing terms yields a solution for hj^{T) 
that is accurate to within 0.1% for all values 0 < / < 
1. Next, as observed in [13, we can use the asymptotic 
expression 


32u{kT) (-1)' 


t Sin(fcr) 
kr 


as kr 00 , 


(62) 


and that C 2 is given by 


C2{k) 


(a - 


1 / 2 ) 


[2/a 


2a 


along with the / = 0 solution hk{T) = hk{T^)jo{kT) to 
infer that the tensor anisotropic stress H^ induces no ad- 
(58) ditional phase shift in hk, so that our earlier expression 


TT 













fraction of free-streaming particles ( f) 


FIG. 1; C 3 is the transfer function factor that accounts for 
the damping of the tensor power spectrum due to tensor 
anisotropic stress. The factor depends on the fraction / of 
the background (critical) energy density contained in free- 
streaming relativistic particles. The figure plots this depen¬ 
dence for 0 < / < 1. 


1S3 for is unchanged. (See da for a complementary 
explanation of this null result, based on causality.) In 
this way, one also sees that 11 ^ damps the tensor power 
spectrum by the asymptotic factor 

C 3 = (63) 

where 

CO 

A=Y,{-ira,^. (64) 

n—0 

For example, keeping the first 4 terms in this sum, we 
find an approximate expression for A: 

(98/3 - 589/2 + 9380/- 55125) 

(15 + 4/)(50 + 4/)(105 + 4/) ’ ^ ^ 

which is accurate to within 1 % for all values 0 < / < 1 . 

If we keep the first 5 terms in the sum, we find an even 
better approximation for A: 

15(14406/^-55770/3-h3152975/2-48118000/-b324135000) 
343(15 -b 4/)(50 -b 4/)(105 -b 4/)(180 -b 4/) 

( 66 ) 

which is accurate to within 0 . 1 % for all values 0 < / < 1 . 
These calculations improve on the accuracy of previous 
calculations mill 

The exact dependence of Cg on / is shown in Fig. ^ 
Note that, as / varies between 0 and 1, the damping 
factor Cg varies between 1.0 and 0.35. In particular, 
if we substitute / = 0.4052, corresponding to 3 stan¬ 
dard neutrino species,the damping factor 0.80313 agrees 
with the results of numerical integrations Hi 113 . When 
the modes probed by the CMB re-enter the horizon, the 
temperature is relatively low (corresponding to atomic- 
physics energies), so we are fairly confident that neu¬ 
trinos are the only free-streaming relativistic particles. 


But when the modes probed by laser interferometers 
re-enter the horizon, the temperature is much higher 
(above the electroweak phase transition, T ^ 10^ GeV), 
so that the physics (and, in particular, the instanta¬ 
neous free-streaming fraction /) is much more uncertain. 
Thus, laser interferometers offer the possibility of learn¬ 
ing about the free-streaming fraction / in the very early 
universe at temperatures between the inflationary and 
electroweak symmetry breaking scales. 

Finally, although Weinberg and subsequent authors 
have concentrated on the tensor anisotropic stress due 
to a single fermionic species (the neutrino), it is straight¬ 
forward to generalize the analysis to include a combina¬ 
tion of species which (*) may each decouple at a different 
time and temperature, and {ii) may be an arbitrary mix¬ 
ture of bosons and fermions. We find that, as long as all 
of these free-streaming species decouple well before the 
modes of interest re-enter the horizon, then all of the re¬ 
sults presented in this section are completely unchanged. 
In other words, in order to determine the behavior of the 
tensor mode function, one only needs to know one num¬ 
ber — the total fraction / of the critical density contained 
in free-streaming particles — even if the particles are a 
mixture of fermionic and bosonic species with different 
temperatures and decoupling times. 


D. Equation-of-state corrections, 

In this section, we consider various physical effects that 
cause the equation of state Wr(z) to deviate from 1/3 dur¬ 
ing the radiation-dominated epoch, and the correspond¬ 
ing modifications that these effects induce in the GWB 
transfer function. Some of these effects have been dis¬ 
cussed previously by Seto and Yokoyama • The devi¬ 
ations 


Swr(z) = Wr(z) — 1/3 (67) 


primarily modify the transfer function through the red- 
shift factor [1 + Zf/) that appears in Ci [see Eqs. lll??ll and 
(ICTll ] ; through the horizon-crossing factor (72 (A) [see Eq. 
(IKHll ]: and through the phase shift [see Eq. ((STJ]. We 
consider here six physical effects which can produce these 
kinds of modifications of the transfer function. 

First, deviations can be caused by mass thresholds in 
the early universe. Suppose that all particle species are 
described by equilibrium distribution functions. Then we 
can write p and p as 


P 

P 


27r2 ^ 
i 


{v? — ^ 2 ) 3/2 
exp[u -y/\±l 
(^2 — 

exp[u-?/J±l 


u^du, 

du, 


(68a) 

(68b) 


where the ith species (with mass m^, and internal de¬ 
grees of freedom) is described by temperature and 
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chemical potential and we have defined the dimen¬ 
sionless quantities = m^/T^ and [^. In 

the denominator, the -I- and — signs are for fermions and 
bosons, respectively. Then the deviation Sw^. is given by 
the exact expression 


Swj. — 5w^ 


(69) 


where 




9*p ' 

represents the contribution from the ith species, 

_ ^ T/ 15 {u^ - ^ 

9^p - / ^ 9i ryn4 4 / f 1 J_ 1 

^ exp[M-j/jJ±l 


(70) 


(71) 


represents the effective number of relativistic degrees of 
freedom, T is conventionally chosen to be the photon 
temperature, and we have defined the function 


f{x,y) 



{v? — 

exp [u — y]±l 


(72) 


species is governed by a classical action that is confor¬ 
mally invariant at the classical level, leading to the usual 
conclusion that the stress-energy tensor is traceless and 
Wr = 1/3. But conformal invariance is broken at the 
quantum level by interactions among the particles, so 
that T/' ^ 0. For example, for a quark-gluon plasma 
governed by SU(IV^) gauge theory, with flavors, and 
gauge-coupling g, the equation of state correction is given 
(up to 0[g^) corrections) by [^1^ 


5w^ 


187r2 (47r)2 2 -h l[N^Nj/{N‘^ - 1)] 


(74) 


Note that this effect can be non-negligible: for large 
gauge groups (be. large in the early universe (prior to 
the electroweak phase transition), the equation of state 
may may easily be reduced from 1/3 by several per¬ 
cent, or more. 

Third, deviations can be produced if the early universe 
behaves like a slightly imperfect fluid. The stress-energy 
tensor for an imperfect fluid contains (in addition to the 
usual perfect-fluid terms) three extra terms whose coeffi¬ 
cients (x, r], and /) represent heat conduction, shear vis¬ 
cosity, and bulk viscosity (see Weinberg [^, Ch. 2.11). 
Of these dissipative effects, only the bulk viscosity term 


For fixed j/^, note that f{x^,yj) vanishes as x^ goes to 
0 or oo; and in between it has a fairly broad peak, 
with a maximum located at and a peak value 

ypeafc _ jjj particular, when yi = 0, then 

the ordered pair {x^^°‘^, is (2.303,1.196) for bosons 
and (2.454,1.125) for fermions. This makes sense: we 

expect 5w^r^ to vanish when Xj <C 1 (since the species is 
relativistic) and when Xj 3> 1 (since the species is non- 
relativistic, and makes a negligible contribution to the 
energy density). In between, when Xj ^ x^ , the fth 
species is cold enough to exhibit non-relativistic behav¬ 
ior, yet hot enough to contribute non-negligibly to the 
energy density. 

Using the above equations, we can compute Sw^(z) 
once we know T^(z) and But let us estimate the 

size of the effect. As a species becomes non-relativistic, 
it produces a maximum equation-of-state deviation 




5fr' 9. Tf 


g.pT^ 


(73) 


in the background equation of state. Furthermore, if Ng 
different species (with the same temperature and similar 
masses) become non-relativistic at the same time, then 
(roughly speaking) the effect is multiplied by Ns (since 
their add). Ultimately, the fractional correction 

Sw^/wj. is model-dependent, but it can conceivably be as 
large as a few percent. 

Second, deviations can be produced by a trace anomaly 
in the early universe. During the radiation-dominated 
epoch, the universe is dominated by highly relativistic 
particles whose masses may be neglected. Thus, each 


+ U^U'')U\^ (75) 

can contribute to the background evolution in an FRW 
universe (see Weinberg |^, Chs. 15.10-15.11). This term 
modifies the continuity equation 


p = -3i7(p + p) + 9Ci7" 


-mp 


l + w — 


SttGCi 
H . 


(76) 


so that, as far as gravitational waves are concerned, the 
effective equation is corrected by 


6wj. 


SttGC 

H 


(77) 


Whereas the three effects discussed thus far produce 
small corrections to the equation of state, it is worth men¬ 
tioning three other effects that can produce much larger 
deviations. The first example is a massive particle species 
that decouples from the thermal plasma before its abun¬ 
dance becomes negligible. Since its energy density falls as 
a~^ (more slowly than the radiation density, which falls 
as a“'*), it can come to dominate the energy density of 
the universe before it decays (if its lifetime Tdecay is suffi¬ 
ciently long). In this case, w drops to zero when the par¬ 
ticle becomes dominant, and rises back to w = 1/3 over 
a timescale given by the decay lifetime Tdecay Second, 
extra-dimensional physics typically modifies the effective 
4-dimensional Friedmann equation. Such modifications 
which, from the standpoint of the GWB, can in some 
cases look like corrections to the effective radiation equa¬ 
tion of state, are primarily constrained by the require¬ 
ment that the Friedmann equation becomes sufficiently 
similar to the ordinary 4-dimensional Friedmann equa¬ 
tion (with ordinary matter) by the time of BBN. Third, 
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due to the weak coupling of the inflaton, the tempera¬ 
ture at the start of the radiation dominated epoch (the 
reheat temperature) can be much lower than the energy 
scale at the end of inflation. In this case, laser interfer¬ 
ometer scales might actually re-enter the horizon during 
the reheating epoch (before the start of radiation dom¬ 
ination), when the equation of state was probably quite 
different from w = 1/3. The actual equation-of-state 
depends on the details of the reheating process, but a 
commonly-discussed value is w = 0, or some value in the 
range 0 < w < 1/3 [^. If w = 0 during reheating, the 
corresponding modification of the GWB might be simi¬ 
lar to the modification due to the long-lived massive relic 
discussed above. 

Note that the first of these six processes can be ex¬ 
pressed as a modification of (/*, and the effects on the 
transfer function can be computed using the methods 
discussed in Ref. M- However, the other five cannot. 


IV. DISCUSSION 

Fig.H illustrates some of the effects discussed in sec¬ 
tion cni In this figure, the solid black curve represents 
the present-day energy spectrum, Ogu,(/, Tq), generated 
by a particular inflationary model — namely, a quadratic 
potential U(</) = (l/2)m^(/^. The red dotted curve il¬ 
lustrates the damping due to tensor anisotropic stress 
from free-streaming neutrinos. We have assumed that 
the free-streaming fraction is / = 0.4052, which is the / 
value for three standard neutrino species which decouple 
around the time of BBN. The green dot-dashed curve rep¬ 
resents the damping due to tensor anisotropic stress from 
various particle species {X particles) which begin free- 
streaming before the scales detected by BBO/DECIGO 
re-enter the horizon and then decay after the scales re¬ 
enter, but prior to the electroweak phase transition. As 
an example, we have assumed that the free-streaming 
fraction is / = 0.5. Finally, the blue dashed curve rep¬ 
resents damping due to a trace anomaly that is present 
above the electroweak scale. For illustration, we have 
assumed that this anomaly, through Eq. m, reduces 
the equation of state from = 1/3 by Sw^. = —0.02. 
This reduction may be achieved by various combinations 
of the number of colors A/, the number of flavors 
and the gauge coupling g; but the point is that we have 
not chosen an unreasonable large value for dw^, given the 
large gauge groups that are often theorized to be present 
at high energies. 

The key point conveyed by Fig. |21is that there are a 
variety of plausible post-inflationary effects that can pro¬ 
duce rather large modifications of the gravitational-wave 
spectrum on laser-interferometer scales, without modi¬ 
fying the spectrum on CMB scales. This is tantalizing, 
since the modifications on laser-interferometer scales re¬ 
flect the primordial dark age between the end of inflation 
and the electroweak phase transition, at energies beyond 
the reach of terrestrial particle accelerators. 



FIG. 2: The black solid curve represents the present-day 
gravitational-wave energy spectrum, Qg^{f,ro), for the in¬ 
flationary model V{(f)) = (l/2)m^0^. The red dotted curve 
shows the damping effect due to (three ordinary mass¬ 
less species of) free-streaming neutrinos. The green dot- 
dashed curve shows the damping effect which arises if free- 
streaming particles make up fifty percent of the background 
energy density at the time rggQ when the modes probed by 
BBO/DECIGO re-enter the horizon. As shown in the figure, 
the particles begin free-streaming sometime before Tggo, and 
decay sometime after TggQ , but prior to electroweak symme¬ 
try breaking. Finally, the blue dashed curve shows the effect 
of a conformal anomaly in the early universe that slightly 
reduces the equation of state from w = 0.33 to w = 0.31 
above the electroweak phase transition. The spectrum will 
also be modified on comoving scales that re-enter the horizon 
during the reheating epoch after inflation; but the range of 
scales affected by reheating is unknown. Finally, note that 
the correlated BBO interferometer proposal claims a sensitiv¬ 
ity that extends beyond the bottom of the figure (down to 
roughly ~ 10“^^) in the frequency range from 10~^ Hz 
to 10° Hz. 


Let us look ahead to future observations. First, as 
measurements of the scalar power spectrum improve, we 
will be restricted to considering inflationary models that 
match precisely the measured scalar perturbation spec¬ 
trum on CMB scales. This is not enough to specify a 
unique model because there is a range of inflaton poten¬ 
tials whose scalar perturbation spectra are indistinguish¬ 
able on CMB scales. For example, two such models are 
shown in Fig.|3 The two curves correspond to two differ¬ 
ent inflationary models whose parameters have been cho¬ 
sen so that they produce the same scalar amplitude and 
tilt on CMB scales, and match the observed scalar ampli¬ 
tude. From the figure, though, we see that their primor¬ 
dial tensor power spectra are distinguishable. The upper 
solid curve corresponds to the simple quartic monomial 
potential V{(j)) = Ai/^, and the lower solid curve corre¬ 
sponds to the axion (or “natural inflation” [^) potential 
V ((/) = Mj [1 -I- cos{(j)/g,)]. 

Note that the two spectra “curve” significantly down¬ 
ward over the broad range of wavenumbers shown in 
Fig.El so that they are not well described by a power-law 
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FIG. 3: The precise predictions of the primordial tensor 
power spectrum for two inflationary models whose parameters 
have been chosen so that they make identical predictions for 
the scalar fluctuation amplitude and spectral index on CMB 
scales. The two models are the quartic monomial potential 
V{(p) = (top), and the axion (or “natural inflation”) po¬ 
tential V{(f)) = M/[1-I-cos(())//i)]. 


(straight-line) approximation over this range. In partic¬ 
ular, a power-law extrapolation based on CMB measure¬ 
ments will give a fairly poor estimate of the real spectrum 
on laser-interferometer scales. In fact, since the real spec¬ 
tra tend to redden increasingly on smaller scales [i.e. as 
the end of inflation approaches), the CMB power-law ex¬ 
trapolation will tend to overestimate the tensor spectrum 
on laser-interferometer scales. 

Second, note that the inflationary consistency relation 
(Pt implies that a gravitational-wave spectrum with a 
larger tensor/scalar ratio r on CMB scales also tends to 
have a “redder” (more negative) tilt nt- For gravitational 
waves with detectably-large amplitudes on CMB scales, 
this leads to the “convergence effect” in Fig.|21 if two dif¬ 
ferent inflationary models produce scalar perturbations 
with the same amplitude and tilt on large (CMB) scales, 
then their tensor spectra tend to approach and even cross 
one another on small (laser-interferometer) scales, as il¬ 
lustrated by the two solid black curves in Fig.|21 As a re¬ 
sult, once one has measured the scalar and tensor pertur¬ 
bations on CMB scales, one expects the tensor spectrum 
on laser-interferometer scales to lie within a relatively 
narrow and predictable range. 

This is both good news and bad news for laser- 
interferometer science. The bad news is that the con¬ 
vergence effect tends to make laser interferometers less 
effective at breaking any remaining degeneracies between 
inflationary models because their predictions for the ten¬ 
sor spectrum are so similar on small scales. The degree to 
which laser interferometers can break such degeneracies 
has been studied in ■ The good news is that the mea¬ 
surements can potentially provide a model-independent 


test of a key inflationary prediction: that the tensor spec¬ 
trum observed at long wavelengths really does extend 
to small wavelengths, as expected from the consistency 
relation Furthermore, precise laser-interferometer 

measurements can serve as sensitive probes of new post¬ 
inflationary physics that modifies the transfer function 
on small scales. 

In this light, let us now imagine some time in the future 
in which we have successfully measured the tensor ampli¬ 
tude on CMB scales, and consider four possible outcomes 
of a laser-interferometer experiment (like BBO). First, 
suppose that BBO detects the gravitational-wave ampli¬ 
tude at the expected value, as exemplified in Fig. 13 This 
would be a quantitative confirmation of inflation gener¬ 
ally, and of the inflationary consistency condition 
in particular, since the large difference in wavenumber 
between CMB and laser-interferometer scales provides a 
large lever-arm to measure the tensor tilt n^. It would 
also suggest an upper bound on the size of any exotic 
effects in the transfer function on laser-interferometer 
scales. 

Second, suppose that BBO detects the gravitational- 
wave amplitude significantly below its expected value. 
Then the interpretation is less straighforward. On the 
one hand, we could interpret the suppression as being 
due to transfer-function effects, such as those discussed 
above. With this interpretation, the suppression would 
imply a rare opportunity to measure physical properties 
of the early universe, at temperatures above the elec- 
troweak scale, when the relevant modes re-entered the 
horizon. On the other hand, we could interpret the sup¬ 
pression to mean that inflation is more complicated than 
we expected: perhaps there is a feature in the inflaton po¬ 
tential which suppresses the primordial tensor spectrum 
on small scales relative to our expectations, or perhaps 
the consistency condition itself is violated (as it is in some 
multi-field models of inflation). Further measurements of 
the tensor spectrum at intermediate scales are needed to 
resolve this issue. 

Third, suppose that BBO detects the gravitational- 
wave amplitude significantly above its expected value. 
At this point it is worth noting that virtually all of the 
transfer-function effects mentioned in section^] (includ¬ 
ing tensor anisotropic stress, the late decay of a massive 
relic species, bulk viscosity, a conformal anomaly, or stan¬ 
dard reheating with w < 1 /3) suppress the gravitational- 
wave spectrum on small scales, relative to large scales. 
In order to enhance the spectrum on laser-interferometer 
scales, one must invoke an even more exotic effect, such as 
a reheating epoch with a low reheat temperature and an 
unusual equation of state [w > 1/3); or the production 
of gravitational waves after inflation due to the collision 
of bubbles from a phase transition; or perhaps extra- 
dimensional physics. 

Fourth, suppose that BBO detects nothing at all. This 
would most likely indicate a fundamental problem with 
inflation, since an period of accelerating expansion pro¬ 
duces a broad and nearly-flat spectrum of primordial 
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tensor perturbations quite generically (regardless of the 
particular model that drives inflation). One way out of 
this conclusion would be to invoke an extreme suppres¬ 
sion effect in the transfer function on BBO scales — e.g. 


due to a reheating epoch with equation of state w = 0 
and a reheat temperature well below 10^ GeV. Determin¬ 
ing which interpretation is correct would require digging 
deeper, with more sensitive experiments on BBO scales. 
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